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Abstract 

We propose a low energy effective field theory of QCD at the scale of pion mass for 
the Nb = 2 sector, Nb being the baryon number, which contains two dibaryon fields 
in addition to the nucleons and pions. It has a well defined counting, is renormalizable 
and the nucleon-nucleon scattering amplitudes are manifestly unitary at leading order. 
We work out a lower energy effective theory for nucleons with energy much lower than 
the pion mass and three momentum comparable to it, which also has a well defined 
counting and is renormalizable. The dibaryon fields must also be kept as explicit 
degrees of freedom in this theory. We calculate the scattering amplitudes at next-to- 
leading order for the 1 5o and 3 Si channels in this framework and obtain an excellent 
description of the phase shifts for center of mass energies in the — 50MeV range. 

PACS: 14.20.Pt, 13.75.Cs, 21.30.Fe, 21.45.Bc, 03.65.Nk . 
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1. Introduction 



Since the original suggestion by Weinberg pQ that the nuclear forces could be understood 
within the framework of effective field theories (EFT) there has been an enormous devel- 
opment of the subject (see [2j [31 HJ E] for recent reviews). A key ingredient of the EFT 
formalism is that the cut-off dependence which is introduced in order to smooth out ultravi- 
olet (UV) singularities can be absorbed by suitable counterterms, and hence any dependence 
on physical scales much higher than the ones of the problem at hand can be encoded in a 
few (unknown) constants. In order to achieve this in a systematic manner counting rules are 
also necessary. 

Weinberg's suggestion consisted of two steps. The first one was calculating the nucleon- 
nucleon (NN) potentials order by order in Chiral Perturbation Theory (%PT) from the 
Heavy Baryon Chiral Lagrangian (HB^L) [6J. The second one introducing the potentials 
thus obtained in a Lippmann-Schwinger (LS) equation. There is no doubt that the first step 
can be carried out within an EFT framework: the renormalized NN potentials are known at 
leading order (LO), next-to-leading order (NLO), next-to-next-to- leading (NNLO) [HE] and 
next-to-next-to- next-to-leading [9], and isospin breaking terms have also been taken care of 
[TO] . These potentials have been evaluated using static propagators for the nucleon fields. 
The use of non-relativistic propagators gives additional contributions starting at two loops, 
the leading order of which have been evaluated in [UJ. The second step however is delicate. 
The potentials obtained in the first step are increasingly singular at short distances as we 
raise the order of they are calculated. Hence the introduction of a regulator in the LS 
equation is compulsory. If a finite cut-off is accepted and allowed to move between a certain 
range, a very successful description of the phase shifts in different partial waves is achieved 
within this approach [HEKE]. Nevertheless, even with the LO potential, it is not clear that 
the scattering amplitude thus obtained is cut-off independent. An alternative to Weinberg's 
approach, which was free from renormalization problems, was proposed by Kaplan, Savage 
and Wise (KSW) [131 El- The key ingredient was to assume that contact interactions are 
enhanced with respect to the standard chiral counting. However, when the NN scattering 
amplitudes were worked out in this approach at NNLO a bad convergence of the series was 
observed, specially in the 3 Si channel [T5] . Since then a number of proposals has been put 
forward [HU [T7] [HI EH EH EEl [221 Ml EH [251 (26] [2ZI [28l [29l [30l EJJ [321 [33] [3H [35l [36l EH 
38j[39, 40j[41j[42j. Let us mention, for instance, ref. [17] in which it is claimed that the 
renormalization program can actually be carried out in the Weinberg approach at LO if the 
potential is expanded about the chiral limit. However, the removal of the cut-off in this 
approach requires unconventional flows and additional counterterms for higher partial waves 
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Here we elaborate on the idea that the difficulties encountered so far in constructing a 
consistent and useful nucleon-nucleon effective field theory (NNEFT) may be a consequence 
of a misidentification of the low energy degrees of freedom. We will assume the NNEFT for 
energy and momentum scales much lower than A x contains two dibaryon fields as explicit 
degrees of freedom, with energy gaps (residual masses) of the order or smaller than the 
pion mass. If the dibaryon fields are naively integrated out, one gets the enhanced contact 
interactions of the KSW approach. We will argue that they must be kept as explicit degrees 
of freedom. 

The relation between dibaryon fields and the KSW approach was noted early |43j . 
Dibaryon fields have also been used in EFT formulations of the three body problem (see 
[4"4"t l4"5l H6] and references there in). However, they have mostly been regarded as a conve- 
nient trick to carry out calculations (see, for instance, [Ej). What it is new in our approach 
is the assumption that they must be included as explicit degrees in the NNEFT □. They 
cannot be integrated out if one wants to keep a natural counting. For a fundamental field 
theory their introduction should be irrelevant, since one can build the appropriate quantum 
numbers of the dibaryon out of the nucleon fields, and their inclusion does not affect the 
symmetries of the theory. For an effective theory, however, where calculations are necessarily 
organized in ratios of scales, it is extremely important to keep the appropriate degrees of 
freedom in the Lagrangian, even if they may appear redundant at first sight. We hope this 
will be illustrate in the paper with sufficient detail. 

We will organize the paper as follows. In the next section we introduce the NNEFT with 
dibaryon fields, and discuss how the calculations must be organized. In section 3 we match 
it to a lower energy effective theory for energies smaller than the pion mass and momenta 
comparable to it. In section 4 we calculate the NN amplitudes at NLO. In section 5 we 
extract the low energy constants from data. In section 6 we critically examine the output 
of the previous section and propose new counting rules. In section 7 we discuss our results 
and section 8 is devoted to the conclusions. 

2. The nucleon-nucleon chiral effective theory with 
dibaryon fields 

We will consider an effective theory for Nb =1 2 sector of QCD for energies much smaller 
than A x about 2m n, being the nucleon mass. The usual degrees of freedom for such a 
theory, namely nucleons and pions, will be augmented by the inclusion of two dibaryon fields, 

1 In a model dependent framework, the inclusion of dibaryon fields as explicit degrees of freedom has 
already been advocated by some authors [1H1 HH1 ESS I5T1I512] ) 
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an isovector {Df) with quantum numbers 1 Sq and an isoscalar (D v ) with quantum numbers 
3 S\. Since ~ A x , a non-relativistic formulation of the nucleon fields is convenient [6j. 
Chiral symmetry, and its breaking due to the quark masses in QCD, constrains the possible 
interactions of the nucleons and dibaryon fields with the pions. The sector without dibaryon 
fields is the standard one PQ, 



C*n =N^ (iD - g A (u ■ |) + ^-)N - ^-{N^Nf - ^{N^Nf+ 



(2.1) 



+ J -f {Tr(d^d"U) + mlTr(U ] + U) } , U = e 2il 



where u 2 = U, = i \ vJ, d^u}, = (d^ + <9 m m]), 7i a is the pion field, r a the isospin 
Pauli matrices, g A ~ 1-25 is the axial vector coupling constant of the nucleon, and f n ~ 
132MeV is the pion decay constant. This is the leading order Lagrangian for the pions 
(0(p 2 )) and the pion-nucleon interactions (0(p)), augmented by the kinetic term of the 
nucleon (which is next-to-leading order for E ~ p ~ m^, but becomes leading order for 
E ~ p 2 /2m N < TOtt). 

The sector with dibaryon fields and no nucleons in the rest frame of the dibaryons reads 

£d = £o(p) + £o( P 2 ), (2-2) 
where C-o{p) is the 0{p) Lagrangian, 

Dt(-ido + 6' ms )D s ] +Dl(-id + 6' mv )D v , (2.3) 



where D s = D a s T a and 5' , i — s,v are the dibaryon residual masses, which must be much 
smaller than A x , otherwise the dibaryon should have been integrated out as the remaining 
resonances have. The covariant derivative for the scalar (isovector) dibaryon field is defined 
as d D s = d D s + d u], D s ]. Co(p 2 ) is the 0(p 2 ) Lagrangian, 

£ 0{p2) = Sl Tr[D s (uM ] u + v)Mv))D\] + s 2 Tr[D\{uM ] u + u*Mu i )D s ]+ 
+ s 3 Tr[DlD s u u ] + s 4 Tr[D s Dlu uo] + ,%Tr[D\D s UiUi\ + 
+ s 6 Tr[D s Dlu iUt } + s 7 Tr[Dlu D s u ] + s s Tr[Dlu i D s u i ]+ (2.4) 
+ VtDl ■ D v Tr[u^Mu ] + uM ] u] + v 2 D\ ■ D v Tr[u u ] + 
+ v 3 Dl ■ D v Tr[uiUi} + v 4 (D^D j + ZW^Tr^], 

M. is the quark mass matrix, which we will take in the isopin limit, namely the average of 
the up and down quark masses m q times the identity matrix. Si, i — 1, .., 8 and Vj, j = 1, .., 4 
are low energy constants (LEC). 
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The dibaryon-nucleon interactions will only be needed at leading order 

C DN =4|(iW 2 rV 2 iV*) J D s , a + ^={N T a 2 r 2 r a N)D\ a + 

V V (2.5) 

+ ^(iWWiV*) ■ D v + ^L(N T r 2 a 2 aN) ■ L>1 
v 2 v 2 

Ai ~ Ax 1/2 , i = s,v. 

The nucleon-nucleon scattering amplitude will be dominated by the dibaryon field. At 
tree level it gives a contribution ~ l/m n A x (for energies ~ m n ) which is parametrically larger 
than the contribution arising from the four nucleon interactions ~ 1/A 2 - The dibaryon field 
propagator gets an important contribution to the self-energy due to the interaction with the 
nucleons (Fig. lb), which is always parametrically larger than the energy E. As a consequence 
the LO expression for the dibaryon field propagator becomes (in dimensional regularization 
(DR) and minimal subtraction (MS) scheme) , 

(2.6) 



I'H TV 



(p = y/m^E) rather than the tree level expression i/(—E + 8' m . — it]). Note that ( 12. 6j) . 
unlike the tree level expression, has the correct positivity properties irrespectively of the sign 
chosen for the time derivative in (12.31) . The unconventional signs for the time derivatives in 
(12. 3p are chosen in this way in order to correctly reproduce the sign of the effective ranges 
later on. They do not imply any violation of unitarity because the correct leading order 
expression for the propagator is H 2 . 6 1) and not the tree level one. From (12. 6ft it follows that 
the leading contribution to the NN scattering amplitude for energies ~ m n is parametrically 
~ 1 / m]/ 2 A^/ 2 , namely slightly suppressed with respect to the tree level estimate for it, 
but still more important than the tree level contribution from the four nucleon interactions 
(~ 1/A 2 ,). The fact that a loop contribution always dominates over a tree level one is kind 
of bizarre in an EFT framework. It probably indicates that the first reliable approximation 
to the true amplitude is the NLO one, namely the first order in which the tree level energy 
dependence is not neglected. It also suggests that NLO contributions (and probably beyond) 
should be resummed in some sort of self-energy. We will see later on that unitarity provides 
the key ingredient to carry out these resummations. 

Equation (12.61) implies that the dibaryon field should not be integrated out unless p <C 5' m . , 
instead of E -C S' m . as the tree level expression suggests. If S' m . <C m^, it should also be kept 
as an explicit degree of freedom in the so called pionless EFT, like in Refs. [531 EH EE]. 

Except for the above mentioned contributions to the self-energy of the dibaryon fields, 
which become LO, the calculation can be organized perturbatively in powers of 1/A X . Hence 
one expects that any UV divergence arising in higher order calculations will be absorbed in 
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a low energy constant of a higher dimensional operator built out of nucleon, dibaryon and 
pion fields (note that the linear divergence in the self-energy of the dibaryon fields due to 
the diagram in Fig. lb can be absorbed in 8 mi ). The renormalized result may be recast in a 
manifestly unitary form [HJ [56]. 

We shall restrict ourselves in the following to energies E ~ m 2 r /A x <C m n , which implies 
nucleon three momenta ~ m n . We shall follow the strategy of [12], which was inspired in the 
formalism of [57] , and shall build a lower energy EFT with no explicit pion fields: the effects 
due to the pions will be encoded in the potentials (and redefinitions of the LECs). We will 
present the calculation of the NN scattering amplitudes at NLO. For these energies the LO 
contribution is ~ l/m 7r A x , which means that we are aiming at including all contributions 
up to 1/AJ. 



3. The lower energy effective theory 

For energies E ~ m 2 r /A x <C m^, the pion fields can be integrated out. This integration 
produces nucleon-nucleon potentials and redefinitions of low energy constants. Since we are 
aiming at a NLO calculation we must keep corrections 0(m„-/A x ) and neglect higher order 
ones. 

In the one nucleon sector pion loops produce corrections which are 0(m^/A x ) and hence 
can be neglected. The same holds true for the dibaryon sectors. However, both single nucleon 
and dibaryon sectors get a contribution 0(m 7r /A x ) from counterterms proportional to the 
quark masses which redefine the nucleon mass and the dibaryon residual masses. 

In the two nucleon sector, the one pion exchange is the only relevant contribution at this 
order, which produces the well known one pion exchange (OPE) potential. 

Pion loops in the dibaryon- nucleon vertices also produce 0(m^./A x ) corrections, except 
for those which reduce to the OPE potential correction to the dibaryon-nucleon vertex, which 
are included in the effective theory and must not be considered in the matching. 

Hence the Lagrangian of the lower energy effective theory at the NLO order reads as 
follows. The sector without dibaryon fields reduces to the LO one in the Weinberg approach 

m, 

d 2 \ C 7 



V Zitln/ 2 2 

+ -iVVW^OK'i^i - x 2 )N^a j rN{x 2 ), 



(3.1) 



where is the OPE potential, 



9a f d 3 Q QiQj ^.(3i-s 2 ) 



V«{x x - x 2 ) = / 2 e ' mX1 ' X2> - ( 3 - 2 ) 

2/| J {27T) 3 q 2 + ml 
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(a) (b) 

Figure 1: a) LO diagram. The shaded circle stands for the resummation of dibaryon selfenergy. b) Dibaryon 
selfenergy diagram. 

The sector with dibaryon fields and no nucleons in the rest frame of the dibaryons reads 

£'d = DU(-ido + 5 ma )D a s + Dl(-id + 5 mv ^)D v , (3.3) 
S mi , i = s, v is the (redefined) dibaryon residual mass, 

5 m a = 8' ms +Am q {si + s 2 ) -25m N 

$m v = 5' mv + irrigV! - 25m N , (3.4) 

SiriN is also proportional to m q and stands for the leading m q corrections to nucleon mass 
(see, for instance, [SB]), which can be reshuffle into 8 m by local field redefinitions. Si, s 2 and 
Vi are the LEC introduced in equation (12.4H . Note that if 6' m . <C m n then the quark mass 
dependence of 5 mi is a leading order effect. The dibaryon- nucleon interactions remain the 
same as in (|2.5p . 

The calculations in this EFT can be organized in ratios E/p and p/A x (recall m n ~ p). 
The UV divergences arising at higher orders will be absorbed by local terms build out of 
nucleon and dibaryon fields. 




4. The scattering amplitudes at NLO 

The scattering amplitudes at LO are given by the diagrams in Fig.l. The corresponding 
amplitudes are, 

1 4A 2 

A\ = - 2 i = s,v, (4.1) 

Subscripts are i = s,v, and superscripts (which appear below) are i = 1 Sq, 3 Si with the 
equivalence s = 1 5 , , v = 3 5*i. We use DR and the MS scheme throughout. These amplitudes 
formally coincide with the leading order of the effective range expansion (ERE) with the 
scattering length given by a, = A ^™ N , and hence they also coincide with the LO amplitudes 



in the KSW approach identifying with i + The bubble resummation chain in 
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(a) 





(b) 



Figure 2: NLO diagrams, a) Diagrams with a one to one correspondence with the ones in the KSW 
approach, b) Extra diagrams that have been considered in our approach. 

the latter plays the role of the dibaryon field. Notice, however, that the use of the dibaryon 
field makes unnecessary the use of the PDS scheme in order to keep a consistent counting 
for p ~ m w . Furthermore, one must keep in mind that 5 mi may contain a leading order 
dependence on the quark masses if 5 mi <C m n , whereas the LO scattering lengths in the 
KSW approach do not depend on the quark masses. 

At NLO we have the diagrams in Fig. 2. The diagrams in Fig.2a have a one to one 
correspondence with the ones in the KSW approach. We have checked that they give the 
same contribution as in ref. [Hj. The diagrams in Fig. 2b are new, and produce the only 
computational difference between the KSW approach and ours at this order. The insertion of 
three momentum dependent four nucleon interaction in the bubble chain in KSW corresponds 
to the tree level energy dependence of the dibaryon field, which becomes NLO in our approach 
as well, and the insertion of a quark mass dependent four nucleon interaction in the bubble 
chain corresponds to the dibaryon mass shift due to (13.41) . which may become LO in our 
approach. Since the diagrams of Fig. 2b do not contribute to the 3 Si — 3 D\ mixing, we 
obtain the same expressions as KSW for it, and hence we will not discuss it further. The 
contributions of Fig. 2b to the scattering amplitudes in the x Sq and 3 Si channels read 



where C s = C$ — 3CV, C v = C$ + Ct, which produce the following extra contributions to 
the phase shifts with respect to the KSW approach 




i = 



(4.2) 




4tt 



I = 



S,V, 



(4.3) 
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5. Low energy constants from data 



In this section we compare the output of our calculation with data and extract the low 
energy constants. This is important for the self-consistency of the approach: if data favors 
5 mi of the order of the pion mass or smaller, then the introduction of the dibaryon field 
makes a lot of sense. If, on the contrary, it delivers 5 mi larger than the pion masses, then 
the introduction of the dibaryon field should be irrelevant. In the first case, if, in addition, 
(5 mj <C m n our approach is expected to produce qualitative differences with respect to KSW, 
at least as far as the quark mass dependence is concerned. 

As discussed in the previous section, our final results for the NLO amplitudes and phase 
shifts may be obtained from those of the KSW approach, by removing the PDS subtractions 
and adding ( 14. 3p . It turns out that this extra contribution ( 14. 2 j) can be written in a form 
which resembles the quark mass insertion in KSW (see (15.21) below). As a consequence, the 
following correspondence between our parameters and, for instance, those of ref. [H] exist 
once the PDS subtractions are removed, 

AA 2 AA 2 1 

^ = C ° ,a = <___ = Cj, (5.1) 

f m , o mi m N o mi 

Note that whereas there is indeed a one-to-one correspondence between our parameters and 
the ones of the KSW approach, the quark mass dependence of these parameters does not 
match, which may become important for eventual extrapolations of lattice data. 

The fits to phase shift data obtained in this way are very unstable, specially for the 1 Sq 
channel. It is convenient to write these amplitudes in a manifestly unitary form [JTJ 156] , 
which is equivalent to the one described above at the order we are working. This is achieved 
by writing, 

A o = T4§72(^-i) 2 i = s,v, (5.2) 
and proceeding analogously for the pion contributions arising from figSi, 

Hi 

\ <5m, I 



n 2 

^ = -^r- 2 ( A -i) i = s,v, (5_ 3 ) 



where Ilj , i — s, v stand for, 



+ -^ln 1 + 4 

4p z V 



(5.4) 



S 



Due to the simple selfenergy structure, it is easy to recast these^contributions into man- 
ifestly unitary expressions, which lead to the following phase shift 



5j = — tan -1 



e , d+Yi, 



S, V. 



(5.5) 



Note that at the order we are working the contribution proportional to Cj can be absorbed 
into a redefinition of 8 mi . This can already be seen at the Lagrangian level. Indeed, a 
local field redefinition of the dibaryon fields of the type D — ► D + cNN, which respects 
the counting (i.e. c ~ A x 3 ^ 2 ), allows to remove the four nucleon terms at the only cost of 
introducing higher order operators [HI 05j H7j. Nevertheless, we will keep these terms for 
the fits in this section because the outcome will be illuminating. The results discussed in 
the following as well as the figures correspond to the formula (I5.5P above. 

5.1 Fitting the phase shifts 

Fitting the expressions (I5.5P to the Nijmegen data for the phase shifts instead of the expres- 
sions in Section 4 leads to much more stable results, and a very good agreement with data at 
NLO. We restrict ourselves to energies in the range — 50MeV, for which our lower energy 
EFT should hold. We present the results of the fit for the LO and the NLO expressions in 
order to keep track of the evolution of the parameters. For the 1 Sq channel at LO we obtain, 

A 2 

= -3.20 ■ 10~ 5 MeV~ 2 , (5.6) 

and the fit is rather poor. At NLO, however, the fit becomes extremely good (Fig. 3), we 
obtain 



C s = -1.225 ■ W- 4 MeV- 2 , A s = 0.0275MeV~ 1/2 , 5 ms = 21.2MeV, 

A 2 (5-7) 
7^ = 3.56 • 10- 5 MeV~ 2 , 

For the 3 Si channel at LO we already obtain a reasonable fit with, 

=2.57- 10- A MeV~ 2 , (5.8) 



2 The unitarization is only approximate for the S\ channel, since coupled channel effects are neglected 

BDI30. 
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Figure 3: The solid line shows the Nijmegen data for the 1 S'o phase shift. The doted curve corresponds 
to fitting at LO the scattering length, while the dashed line corresponds to fitting at NLO the scattering 
length, the effective range and the first shape parameter (v2 S )- The dot-dashed curve corresponds to fitting 
S in the energy range shown. 
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which becomes better at NLO (Fig.4) with, 



C v = -0.228 • lO^MeV' 2 , A v = 0.0332Mel/" 1/2 , 5 mv = 32.0MeV, 

A 2 (5-9) 
— v - = 0.344- W- 4 MeV- 2 } 

We observe that the values for 5 mi are small, namely 5 mi ~ m 2 /A x <C m n . This explains 
the large variations of A 2 /5 mi in going from LO to NLO: 5 mi — 5' m . is not a small correction 
to 5' m . in (13.41) but a quantity of similar size. A 2 take reasonable values ~ VAx ~ 
and C5 and CV take values larger than expected with the size assignment ~ 1/A 2 . In 
fact, a size assignment ~ 1/A x m n appears to be more appropriate. Note that such a size 
assignment would not allow to remove these terms by a local field redefinition which respects 
the counting, like the one discussed after (15. 5p . In fact, if one sets C s = 0, the good fit to 
data for the 1 S'o channel is spoiled. We will return to this point in section [671 

If we wish to analyze the convergence of the EFT results for the phase shifts, fitting the 
LECs to data at each order may not be the optimal way to proceed. Since the EFT should 
work the better the lower the energy is, it appears reasonable to us to extract the LECs from 
the ERE, as advocated by some authors [24, 25J. 



5.2 Inputting the ERE parameters 

At LO the scattering lengths a s = — 23.7 fm and a v = 5.42/m are sufficient as an input. We 
obtain 

A 2 



-1.61 • 10~ a MeV- 



J m s 



A 2 

v -= 3.68 • lO^MeV- 2 . 



(5.10) 



At NLO we need in addition the effective ranges r s = 2.67 fm and r v = 1.83/m and the 
shape parameters V2 S = — 0.476/m 3 and v^v = — 0.131/m 3 [59j. We obtain, 

C s = -0.874 • 10" 4 Me1/- 2 , A s = 0.0239MeV^ 1/2 , 5 ms = 18.9MeV 

A 2 (5-11) 
— ^ = 3.02 • lO- 5 MeV~ 2 , 

Orn 3 

C v = -0.768 • 10- 4 MeV~ 2 , A v = 0.0292MeV~ 1/2 , S mv = 19.0MeV 

A 2 (5-12) 
" = 4.49 ■ 10- 5 MeV~ 2 . 



From the plots in Figj3] and FigJH we see that inputting the parameters from the ERE 
produces a less satisfactory description of data than the fits. Nevertheless, it might provide 
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Figure 4: The solid line shows the Nijmegen data for the 3 Si phase shift. The doted curve corresponds 
to fitting at LO the scattering length, while the dashed line corresponds to fitting at NLO the scattering 
length, the effective range and the first shape parameter (v2t)- The dot-dashed curve corresponds to fitting 
S in the energy range shown. 

more convenient procedure to analyze the convergence of the series. It also produces smaller 
values of 5 mi , which make the contributions of the dibaryon fields more important at low 
energies. 

6. Rethinking the counting 

The numbers for 5 mi , i — s,v, of tens of MeV, clearly indicate that the dibaryon fields must 
be kept as explicit degrees of freedom essentially at all energies (only for E <C 8 2 /A x ~ 
lMeV it is justified to integrate them out). In addition, since 5 mi ~ m 2 /A x , it implies that 
8' m . < m^/A x , and hence of size comparable to the NLO contributions. This suggests that 
the actual LO of the dibaryon field propagator should be ir/A 2 m^p for p ~ m n rather than 
(12. 6p . This expansion will fail at very low energies, but the manifestly unitary version of it 
(beyond leading order) is expected to produce sensible results in the very low energy region 
as well. This can be checked from our expressions (15. 5ft by keeping only the dependences in 
5 mi which gives rise to LO corrections in the above counting (this is dropping the terms with 
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Ci and keeping only the first term in (I5.4p ). namely, 



5i 



tan 



A^m N p 



(6.1; 



Fits to data in this case, shown in Figj5j are only slightly worse than the ones displayed in 
Figj3] and FigJH They deliver, 



.4, 



0.0305Me\/- 1/2 
0.0365MeV- 1/2 



-WAMeV 



OAOOMeV 



A 2 

— i = -4.82 • 10- 5 MeV~ 2 



(6.2) 



0.00333MeK" 



From the numbers above, it is clear that this expansion is expected to work better for the 
3 Si channel than for the 1 5'o one. It is also interesting to note that the numbers obtained 
for Ci from the fits in sectional namely Cj ~ l/A^m^, produce in (15. 5p the expected size 
for a correction 0(m n / A x ) to the formula (16.11) above. Then the Ci in (15.51) appear to be 
simulating next order corrections to the phase shifts. This is also so for the second term in 
(15. 4p (the last terms in it give rise to a next-to-next order corrections). For the 1 Sq channel 
this is specially important: if one sets C s = in ( 15. 5ft but retains the second term (or the full 
expression) in (15.41) . no good fit to data is achieved. Hence, both terms must be consistently 
neglected (or taken into account). 

In order to evaluate the impact of recasting the amplitude in a manifestly unitary expres- 
sion, let us expand the formula (16.11) around 5 mi = up to NLO terms, i.e., tan _1 (x) ~ | — -, 
and fit this expression to data in the energy range 2 — 50MeV. The expanded expressions 
still produce very good fits in this energy range (see Fig. [5]) and deliver similar values for 
the parameters, 



.4, 



0.0343Me\/- 1/2 
Q.0367MeV- 1/2 



-23.8MeV 



OAOAMeV 



A 2 

= -4.96 • 10- 5 MeV~ 2 



A 

S j n 



(6.3) 



0.00333MeV 



As mentioned before, this expansion breaks down at very low energies. The manifestly 
unitary formula (16. lh essentially fixes up this failure. Alternatively, at these low energies 
one should change the counting to p ~ 5 mi m^. Then one can match our lower energy 
EFT to an even lower energy EFT with only contact interactions between nucleons (rather 
than finite range potentials like Yukawa's), and between nucleons and dibaryons. In fact 
it is within this very low energy EFT that dibaryon fields have mostly been considered in 
the literature [1U H7J It would be interesting to have the explicit relation between the 
parameters in our fundamental theory and the ones in this very low energy EFT, which is 
left for future work. 
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Figure 5: The solid line shows the Nijmegen data for the phase shifts 1 5o (Fig. a) and 3 5i (Fig.b). The 
dashed and doted curve correspond to the fit in the range 2 — 50MeV of (|6.1[) and its expanded version (as 
explained in the text) respectively. The expanded expression fails to reproduce the data at very low energies 
while fits better than the full expression at p ~ m,. Notice that the expanded expressions for the 1 5'o( 3 5'i) 
channel tend to — oo(oo) as the energy tends to zero. 
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7. Discussion 



Let us first elaborate on the apparent contradiction between our hypothesis that dibaryon 
fields must be considered fundamental degrees of freedom and Weinberg's claim that the 
deuteron is not an elementary particle [60]. The point is that one should not identify the 
3 S\ dibaryon with the deuteron. We can substantiate this statement by computing the value 
of Z — | (Q\ D v \d) | 2 (d stands for the deuteron), the projection of the dibaryon field on the 
deuteron, from the LO dibaryon propagator: 

/ d*xe^ (0| T{Di(0)Di\x)} \Q) \ E=p0 ^ Ed = , + . . . , (7.1) 

where Ed is the binding energy of the deuteron. At LO 

/+ ifiv 
d A xe ip - x (ft | T{Dl(0)Di\x)} \Q) = T% = Ti . (7.2) 

from which one easily obtains 



r A 2 m Ny /-m N (E+ie) . , 

T °m~\ ~ O m + \O m \ (7 Q\ 

hm E^E d j^r = A4mlr « 1. (7-3) 



Hence the projection of the dibaryon field on the physical state, i.e. the deuteron, is para- 
metrically small, and hence there is no contradiction between Weinberg's statement that the 
deuteron is mainly a nucleon-nucleon bound state, and considering dibaryons as basic de- 
grees of freedom. This is possible because the interaction of the dibaryon with the nucleons 
is a leading order effect. 

Let us next discuss the quark mass dependence of the scattering lengths a in our approach. 
At leading non-vanishing order they read 

— — ~ c + c'm q + c"m q In m q , (7.4) 

(c, d and c" stand for quantities which do not depend on the quark masses) in which all 
three terms are equally important (the chiral log arises from the counterterm of the rightmost 
diagram in Fig. 2a). The most important correction to this formula comes from a particular 
momentum region, uncovered in [61] and recently discussed in [11] , of two loop self-energy 
diagrams of the dibaryon fields which produces terms proportional to m^ 4 . Hence, the 
lattice results for the scattering lengths |62j [63] are expected to be very sensitive to the 
quark masses used. 

Let us finally mention that the inclusion of dibaryon fields as fundamental degrees of 
freedom may help understanding certain enhancements in iV-body forces. In the Weinberg's 



15 



original counting [T] A-body forces were suppressed by powers of A^~ 3N . However, the 
introduction of dibaryon field allows to introduce local terms in the baryon number sector N 
which are only suppressed by A x 3N ^ 2 if N is even or A x 5 3N ^ 2 if ]\f [ s odd, which indicates 
that A-body forces are expected to be enhanced with respect to Weinberg's counting. In 
order to illustrate it, consider, for instance, a term rD^N^ND, r ~ 1/A^- If the dibaryon 
field is integrated out and the outcome is expanded considering the momenta small, a six 
nucleon contact term r'N^N^N^NNN, r' ~ 1/(A X 5^J is induced, which is indeed enhanced 
with respect to the Weinberg's counting (5 m <C A x ). It is interesting to note that this term 
is of the same 'unnatural' size as the one found in [54], if the cut-off A and the coupling 
constant g are taken of natural size, namely A ~ A x and g 2 ~ 1/A X , in that reference. 



8. Conclusions 

We have proposed a NNEFT at the energy scale of the pion mass, which has two dibaryon 
fields as explicit degrees of freedom in addition to the nucleons and pions. We have matched 
it to a lower energy effective theory in which the pion fields have been integrated out. We 
have pointed out that different counting rules are required for p < m 2 T /A x and for p ~ m n , 
and have focused on the latter case. Both effective theories are renormalizable and give rise 
to manifestly unitary amplitudes at leading order. We have calculated these amplitudes at 
NLO and showed that they produce a very good description of the phase shift in the energy 
range 5 — 50MeV. Once unitarized the good description extends to the full range — 50MeV. 
The residual masses delivered by the fits are small, which indicates that the dibaryon fields 
must also be kept as explicit degrees of freedom at low energies. 
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